Comparative study of electrical characteristics for n-type 4H–SiC planar and trench MOS capacitors annealed in ambient NO
Shen Zhan-Wei1, Zhang Feng1, 3, †, Dimitrijev Sima2, Han Ji-Sheng2, Yan Guo-Guo1, Wen Zheng-Xin1, Zhao Wan-Shun1, Wang Lei1, Liu Xing-Fang1, Sun Guo-Sheng1, 3, Zeng Yi-Ping1, 3
Key Laboratory of Semiconductor Material Sciences, Institute of Semiconductors, University of Chinese Academy of Sciences, Beijing 100083, China
Queensland Micro- and Nano-technology Center, Griffith University, Nathan 4111, Australia
College of Materials Science and Opto-Electronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China

 

† Corresponding author. E-mail: fzhang@semi.ac.cn

Project supported by the National Basic Research Program of China (Grant No. 2015CB759600), the National Natural Science Foundation of China (Grant Nos. 61474113 and 61574140), the Beijing NOVA Program, China (Grant No. Z1611000049161132016071), China Academy of Engineering Physics (CAEP) Microsystem and THz Science and Technology Foundation, China (Grant No. CAEPMT201502), the Beijing Municipal Science and Technology Commission Project, China (Grant Nos. Z161100002116018 and D16110300430000), and the Youth Innovation Promotion Association of Chinese Academy of Sciences (Grant No. 2012098)

Abstract

The interface properties and electrical characteristics of the n-type 4H–SiC planar and trench metal–oxide–semiconductor (MOS) capacitors are investigated by measuring the capacitance voltage and current voltage. The flat-band voltage and interface state density are evaluated by the quasi-static method. It is not effective on further improving the interface properties annealing at 1250°C in NO ambient for above 1 h due to the increasing interface shallow and fast states. These shallow states reduce the effective positive fixed charge density in the oxide. For the vertical MOS capacitors on the () and () faces, the interface state density can be reduced by approximately one order of magnitude, in comparison to the result of the planar MOS capacitors on the (0001) face under the same NO annealing condition. In addition, it is found that Fowler–Nordheim tunneling current occurs at an oxide electric field of 7 MV/cm for the planar MOS device. However, Poole–Frenkel conduction current occurs at a lower electric field of 4 MV/cm for the trench MOS capacitor. This is due to the local field crowded at the trench corner severely causing the electrons to be early captured at or emitted from the SiO2/SiC interface. These results provide a reference for an in-depth understanding of the mobility-limiting factors and long term reliability of the trench and planar SiO2/SiC interfaces.

1. Introduction

Silicon carbide is a promising candidate for advanced power switching devices, owing to the superior material properties such as high thermal conductivity, high critical electric field, and the larger carrier saturation velocity.[14] The full potential of SiC metal–oxide–semiconductor field-effect transistor (MOSFET) is still far from being reached because of the relatively high interface state density (Dit) at the SiO2/SiC interface.[510] The improvement of the channel mobility can be achieved via the post-oxidation annealing (POA) in NH3,[11] N2O,[12] Ar, POCl3,[8, 13] and sodium[14] atmosphere. Among them, nitridation by nitric oxide (NO) is considered to be an effective industrial technique to improve the interface properties. The interface state density can be reduced to a low value for the MOS capacitors annealed in NO at 1250 °C. Nonetheless, the conventional high-frequency method cannot detect the fast states responding to higher than 100 MHz. It is believed that NO annealing at a temperature above 1250 °C is not recommended due to the generation of very fast states.[15, 16] The variation of fast interface state with NO annealing time at 1250 °C or below has not been investigated in detail. Meanwhile, the effects of fast and shallow states on the flat-band voltage shift and channel mobility of the MOSFETs are still ambiguous.

On the other hand, fundamental studies of SiC planar MOS capacitors formed on () and () faces have been performed in the past few years. It is proved that MOSFETs fabricated on these nonpolar faces with NO annealing exhibit higher channel mobilities than those fabricated on (0001) face.[1719] Nevertheless, the channels are formed on a sidewall of the trench for the trench MOSFET. The capacitance–voltage measurement is therefore difficult to use for directly evaluating the interface quality of the sidewall (we call it vertical MOS capacitor in this paper). Anisotropic channel properties on the trench sidewall were mostly evaluated through the value of field effective mobility of the fabricated trench MOSFET.[1921] In addition, the most commonly cited intrinsic oxide degradation mechanism is the Fowler–Nordheim tunneling for the SiC planar MOSFETs.[22, 23] Different mechanisms of charge transports at the SiO2/SiC interface have not been well understood. Therefore, comparisons of electrical characteristic and interface property between the trench and planar MOS capacitors are required in order to obtain a comprehensive understanding of the anisotropic channel mobility and reliability issues for the 4H–SiC MOSFETs.

In this study, the interface properties of planar MOS capacitors for different nitridation times are presented. We also propose a simple method to estimate the density of interface states of the vertical MOS capacitor by using the quasi-static capacitance voltage measurement. Comparison of electrical characteristic between the planar MOS and trench MOS capacitor is discussed in detail.

2. Experiments

Figure 1(a) shows the schematics of one unit cell 4H–SiC planar and trench MOS capacitors, which were fabricated on the n-type 4° off-angle silicon face (0001) epitaxial substrates. The wafers were purchased from TYSiC, with a bulk doping concentration of approximately 8 × 1015 cm−3. For samples A, B, and C, the wafers were prepared for planar MOS capacitors. For sample D, 2-μm deep and 3-μm wide trenches were obtained by inductively coupled plasma (ICP) etching using mixed gases of SF6, O2, HBr. Stripe-shaped unit cell trenches were patterned on the wafer with the SiO2 mask. The directions of the trenches are perpendicular and parallel to the primary flat of the wafer to form the vertical MOS capacitors on the () and () faces respectively, as shown in Figs. 1(a) and 1(b). All the samples were dipped in buffered HF (BOE) after solvent and standard RCA cleaning process. Subsequently, samples A, B, and C were subjected to dry oxidation at 1250 °C for 2.5 h. Post-oxidation annealing (POA) was subsequently performed in NO ambient at 1250 °C for 1 h, 2.5 h, and 5 h, separately. These samples are marked as NO_xh as listed in Table 1 (where x means that the sample is annealed in NO for x hours after dry oxidation). In addition, sample D was annealed in NO at 1250 °C for 2.5 h after dry oxidation of 75 min. A flat reference capacitor was also fabricated on the substrate of sample D. Aluminum of 300-nm-thick was subsequently evaporated and patterned using standard photolithographic procedures to form the gate electrode for each of the samples. BOE was used to remove the oxide grown on the back side, followed by the formation of back contact using nickel.

Fig. 1. (color online) (a) Schematics of one unit cell trench and planar MOS capacitors and (b) definitions of the wafer plane orientation and the trench MOS direction.
Table 1.

Oxide thickness of all the samples obtained by accumulation capacitance (CV) and SE (trench sidewall by SEM).

.

The high frequency (10 kHz) capacitance voltage (HFCV) and the quasi-static capacitance voltage (QSCV) measurements were done using Agilent B1500A. Some useful information can be obtained by the DC and HFCV measurement to set the various QSCV parameters such as the integration time, current range, leakage current and the capacitance value appropriately. For most of the cases, the integration time needs to be longer than the time constant of the QSCV response. To obtain an accurate value of the quasi-static capacitance, it may be necessary to adjust the parameters and measure them by trial and error through comparing the data. The QSCV measurement was performed in steps of 0.2 V for the gate voltage sweep and in an integration time of 1 s for the charge calculation. The interface state density was extracted by the proposed method i.e., the “CψS method”.[15] Meanwhile, we used the “Terman method”[24] to evaluate the interface state density in the planar MOS capacitor. Current–voltage (IV) measurements have also been performed to confirm the process quality and the extrinsic failures of oxide. All measurements were done at room temperature.

3. Results and discussion

Table 1 summarizes the oxide thickness determined by CV accumulation capacitance. The tOX of the vertical MOS capacitor was achieved by using scanning electron microscopy (SEM). These values are in good agreement with the results obtained by spectroscopic ellipsometry (SE) for all the samples. The thickness of SiO2 increases slightly with the increasing nitridation time as can be seen by comparing samples A, B, and C. However, due to the anisotropic behavior of oxidation, the tOX of the trench sidewall is almost three times larger than that of the horizontal face for sample D. The tOX determined by SE for all the samples is used to extract the flat-band voltage (VFB), effective oxide charge density (Qeff), and the interface states density (Dit) as discussed below.

3.1. Planar MOS devices

The normalized CV curves for the planar MOS devices are extracted from QSCV and HFCV measurements by sweeping the gate voltage from depletion to accumulation. The flat-band voltage (VFB) shift as shown in Fig. 2(a), results from the effective oxide charge density (Qeff) at the SiO2/SiC interface. The Qeff is the sum of the interface states density (Dit), fixed oxide charge density, oxide trapped charge density and mobile charges density.[25, 26] The expression for the Qeff can be written as

where ms is the work function difference between SiC (3.81 eV) and Al metal (4 3 eV), COX is the oxide capacitance per unit area. Listed in Table 2 are the VFB and Qeff obtained by the two methods. The VFB obtained by the QSCV method is left shifted in comparison to that obtained by the HFCV method as shown in Fig. 2(a). Furthermore, the curves of semiconductor capacitance (CD) versus gate voltage for sample C obtained by the two methods are shown in Fig. 2(b). It is found that there is a large difference between CD and the theoretical semiconductor capacitance CD, theory at a given gate voltage by using the QSCV method. Therefore, the VFB is underestimated by the conventional HFCV method, because the fast and shallow states also contribute to interface trap capacitance (Cit).[16, 27] These interface trap capacitances are not measured accurately, leading to a wrong value of the flat capacitance.

Fig. 2. (color online) CV curves of (a) planar MOS capacitors for different NO annealing times, and (b) semiconductor capacitance CD versus gate voltage for sample NO_2.5 h obtained by QSCV and HFCV method. The CD, theory versus gate voltage is also plotted.
Table 2.

Flat-band voltages and effective oxide charge densities determined by the QSCV and HFCV method for the planar MOS capacitors for different NO annealing times.

.

The values of Qeff determined by the QSCV method (Qeff,QSCV) are 2.35 × 1011, 2.82 × 1011, and 1.30 × 1011 cm−2 respectively for the planar MOS capacitors annealed in NO atmosphere for 1 h, 2.5 h, and 5 h. In the previous reports, it was demonstrated that the NO annealed oxides possess N-induced hole traps, which are positively trapped charges and scaled with the NO annealing time.[28, 29] Nevertheless, in this experiment, Qeff,QSCV for the sample NO_5 h is smaller than that for the sample NO_2.5 h, and Qeff,HFCV even switches from positive to negative values. We speculate that the shallow state density should increase with the increase of NO annealing time. These shallow states (donor-type traps) exist near the SiO2/SiC interface in the oxide and trap more electrons, leading to a decrease of effective positive fixed charge density. The density of shallow states probably increases with a slow rate when NO annealing time ranges from 1 h to 2.5 h. Thus, it has almost no effect on the total Qeff. A sharp increase in the density of shallow states occurs when NO annealing time ranges from 2.5 h to 5 h, which results in reductions of Qeff,QSCV and Qeff,HFCV. As shown in the last row of Table 2, the values of effective interface fast state density (Qit,eff) are estimated by subtracting Qeff,HFCV from Qeff,QSCV. Note that Qit,eff increases in proportion to the NO annealing time, which may be related to very fast states existing near the SiO2/SiC interface in SiC (the fast states will be discussed later). Therefore, the channel mobility is still below 50 cm2/V·s, which is probably caused by a reduction of the free carrier density at the SiO2/SiC interface and by the scattering of the trapped electrons.[30, 31]

The interface state density is evaluated by using the QSCV (Dit, QSCV) and HFCV (Dit, HFCV) method. The former method was presented in detail in Section 3, together with being used for calculating the vertical MOS capacitors. Figures 3(a)3(d) show Dit,QSCV, Dit,HFCV, and Dit, > 10 K for the planar MOS capacitors annealed in NO for 1 h, 2.5 h and 5 h, respectively. Here, Dit, > 10 K is for the fast state capacitor responding to higher than 10 kHz, and is achieved by subtracting Dit,HFCV from Dit,QSCV. Sample NO_1 h possesses the smallest Dit,QSCV as shown in Fig. 3(a). Meanwhile, the interface state density is beyond a value of 1013 cm2·eV−1 near the conduction band edge for 5 h NO annealing time. As shown in Fig. 3(b), there is no obvious difference in Dit, HFCV with increasing NO annealing time due to the underestimation of the fast state. It is noted that Dit,>10 K does not scale linearly with nitrogen incorporation as indicated in Figs. 3(c) and 3(d). Figure 3(d) also shows that Dit,10 K at ECET = 0.3 eV is reduced by approximately 86% for the sample annealed in NO for less than 2.5 hours.

Fig. 3. (color online) Distributions of (a) Dit, QSCV, (b) Dit, HFCV, (c) Dit, > 10 K and (d) fast states density at ECET = 0.3 eV.

Based on the above-mentioned evaluation results, it is concluded that NO annealing at 1250 °C for above 1 h is not effective in further improving the interface properties due to the increasing fast and shallow state density at the SiO2/SiC interface. These fast states can respond to the frequencies higher than the measurement frequency of 10 kHz in this study and are almost saturated for a longer NO annealing time. However, the interfacial electron trapping in the shallow state will occur at the SiO2/SiC interface. These shallow states are donor-like, thus being neutralized after capturing the electrons. The density of shallow states is rapidly increased when NO annealing time is increased beyond 2.5 h, which causes a reduction of Qeff. The fast and shallow states both determine the flat-band voltage of the MOS capacitor and restrict the improvement of the channel mobility of SiC MOSFET.

3.2. Trench MOS devices

We investigate the crystallographic orientation dependence of the interface state density of the trench MOS capacitor annealed in NO at 1250 °C, especially for () and () faces devices. An appropriate thickness of SiO2 serves as the dry etching mask, which is annealed in O2 at 1000 °C for the densification. Therefore, the surface roughness of the trench sidewall may be reduced. Subsequently, the etched substrate is thermally oxidized to form a sacrificial oxide with a thickness of approximately 20 nm. The sacrificial oxide layer is removed by diluted HF. Thus, the defect density on the SiC surface introduced by the previous etching process may be reduced. Some parts of the sample are patterned by photoresist, to produce a planar MOS capacitor at the same substrate without any vertical edge as shown in Fig. 1(a). The quasi-static capacitance for the trench MOS structure (CQS,trench) is given by the sum of three MOS capacitors as follows:

where C2 is the vertical MOS quasi-static capacitance per unit area for the () or (); C1 and C1′ are the planar MOS quasi-static capacitances per unit area for the (0001) faces; S1, S1′, and S2 are the areas of gate electrode for the three MOS capacitors respectively. Here, C1 can be approximated as C1′, because the CV curves for the two capacitors show a similar appearance.

Based on the quasi-static equivalent circuit shown in Fig. 4(a), the relationship between the quasi-static capacitances and interface trap capacitances is given by:

Combining Eqs. (2) and (4) yields
where CD − 1, COX − 1, and Cit − 1 are the semiconductor, oxide and interface trap capacitances per unit area for the planar MOS capacitor, respectively. Likewise, CD − 2, COX − 2, and Cit − 2 are those for the vertical MOS capacitors, respectively. The oxide capacitance of the vertical MOS capacitor can be determined by the accumulation, which is based on the equivalent circuit shown in Fig. 4(b). The interface state densities for the planar and vertical MOS capacitors can be calculated from
The flat-band voltage of the fabricated MOS capacitor is roughly estimated from the surface potential (ψS) and gate voltage curve as shown in Fig. 5(a). For the vertical MOS capacitors fabricated on the () and () faces, the values of VFB are −4.4 V and −4.2 V, and the values of Qeff calculated from Eq. (1) are 6.85 × 1011 and 6.60 × 1011 cm−2, respectively. As we can see, the vertical MOS capacitor has a large negative flat-band voltage shift. This is because the surface roughness of the trench sidewall causes the generation of positive fixed charge after thermal oxidation and NO annealing. Meanwhile, due to the metal ion contamination of the substrate surface in the ICP process, the calculated Qeff of 7.33 × 1011 cm−2 is larger than the previous results for the planar MOS capacitors under the same NO annealing condition.

Fig. 4. (color online) Equivalent circuits for the trench MOS capacitor in (a) depletion using QSCV method and (b) in strong accumulation using HFCV method, where C1, CD − 1, COX − 1, and Cit − 1 are the quasi-static, semiconductor, oxide and interface trap capacitances per unit area for the planar MOS capacitors, respectively. C2, CD − 2, COX − 2, and Cit − 2 are those for the vertical MOS capacitors, respectively.
Fig. 5. (color online) (a) Comparison among the curves of surface potential versus gate voltage curves of 4H–SiC MOS capacitors on (), (), and (0001) faces, and (b) distributions of Dit extracted from high-frequency QSCV method for those MOS capacitors, with Dit,11 − 20, Dit,1 − 100, and Dit,0001 corresponding to the MOS capacitors on (), (), and (0001) faces, respectively.

Figure 5(b) shows the interface state densities , , and Dit,0001 for the MOS capacitors on (), (), and (0001) faces on the same substrate, respectively. The slight difference between and is caused by the locally nonuniform oxide thickness due to the oxidation anisotropy,[32] which further increases the interface state density at the SiO2/SiC interface. It has been found that the value of Dit at the edge of the conduction band increases as a function of oxide thickness due to the existence of deep acceptor-type surface states.[33] Even though the tOX of the vertical MOS capacitor is larger than that of the planar MOS capacitor, nitridation is carried out under the same NO condition after oxidizing both of the MOS devices. By using the aforesaid method of calculation, it is directly demonstrated that and are smaller than Dit,0001 due to the more superior MOS interface quality of the () and () faces.[30] In spite of the thicker oxide, the interface trap density for the () face at ECET = 0.3 eV can be reduced by approximately one order of magnitude in comparison to that for the (0001) face.

3.3. Comparisons of current–voltage characteristic between planar and trench MOS capacitors

Figure 6(a) shows FN plots for the fabricated planar MOS capacitor, which can be expressed as follows:[34]

where me is the free electron mass, q is the electronic charge, h is the Planck’s constant, and ϕB is the barrier height between oxide and semiconductor. Using the slope of the curve of ln(I/E2) versus 1/E as shown in Fig. 6(a), ϕB can be calculated with a known electron effective mass mox in oxide (0.5me). The barrier height is 2.7 eV, which is close to the value of the band offset between 4H–SiC and SiO2 in a previous study.[35]

Fig. 6. (color online) (a) Electric-field-dependent Fowler–Nordheim (FN) IE characteristics and (b) Poole–Frenkel (PF) IE characteristics for the planar MOS and trench MOS capacitors respectively, with insets showing the energy-band diagrams for electron emission in both of the SiC MOS capacitors. Simulated two-dimensional (2D) electric field distribution when oxide breakdown occurs for the trench MOS capacitors is plotted.

In contrast, a different current–voltage characteristic is obtained for the trench MOS capacitor. A linear function between ln(I/E) and E1/2 is achieved as shown in Fig. 6(b). The current density associated with Poole–Frenkel emission is determined by[36, 37]

where εr is the relative dielectric permittivity, k is the Boltzmann’s constant, ε0 is the permittivity of free space, T is the temperature, q is the electronic charge, Nc is the conduction-band effective density of states, μ is the electron mobility, and ϕt is the potential barrier height for electron emission from the trapped state in oxide. In Eq. (10), the slope of ln(J/E) versus E1/2 can be written as
The relative dielectric permittivity of SiO2 can be obtained from Eq. (11) as follows:
The value of kslope is calculated to be 0.01448 from a linear fitting curve as shown in Fig. 6(b). Thus, the relative dielectric permittivity εr = 4.0 can be extracted from Eq. (12), which is in good agreement with the value for SiO2.

As indicated in the energy-band diagram in the insert of Fig. 6(b), the field emission of electrons will dominate the current transport for the trench MOS capacitor. This is due to the trapped positive charge existing in the thick oxide corner of the trench. Furthermore, as shown by the simulation result of 2D electric field distribution in Fig. 6(b), the local field crowding will severely cause the local electrons to be captured in or emitted from the interface traps. Thus, the PF conduction current can take place at a critical electric field of EC = 4 MV/cm for the trench MOS device. This value is lower than that of planar MOS device (FN tunneling current occurs at EC = 7 MV/cm), probably because more uniform thickness of oxide is formed in the planar MOS capacitor. It is thus concluded that the peak value of the electric field in the gate oxide of trench MOSFET should be reduced to below 4 MV/cm, for the long term reliability and threshold voltage stability.

4. Conclusions

The capacitance voltage and current voltage are measured to evaluate the flat-band voltages, effective oxide charge densities, the interface state densities and oxide breakdown characteristics of the trench and planar MOS capacitors. We report for the first time in this paper the achievement of Dit by measuring the quasi-static capacitance voltage, for the trench MOS capacitor in which only two mask layers are used.

The flat-band voltage is severely underestimated by the conventional HFCV method. Meanwhile, the interface state density does not decrease significantly with the NO annealing at 1250 °C for above 1 h, due to the increasing shallow and fast states generated by the nitrogen incorporation. These shallow states result in a decrease of effective positive fixed charge density because more electrons are trapped into the states near the SiO2/SiC interface in the oxide. On the other hand, it is found that Dit at ECET = 0.3 eV for vertical MOS capacitor on the () face can be reduced to below 1012 cm2·eV−1, which is smaller than that for planar MOS capacitors on the (0001) face. However, the electric field in the gate oxide of trench MOS capacitor should be reduced to below 4MV/cm due to local field crowded at the corner of the trench gate. It is concluded that the trapped positive charges at the trench oxide result in the electrons captured or emitted early. Therefore the PF conduction current occurs at the lower electric field, in comparison to the FN tunneling current for the planar MOS capacitor. The present results will provide a better understanding of the interface properties and reliability of oxynitrided 4H–SiC MOSFETs.

Reference
[1] Cooper J A Melloch M R Singh R Agarwal A Palmour J W 2002 IEEE Trans. Electron. Dev. 49 658
[2] Zhang F Liu S B Dong L Liu X F Fan Z C Liu B Yan G G Wang L Zhao W S Sun G S He Z Yang F H 2013 Chin. Phys. 22 097302
[3] Chen S Z Sheng K 2014 Chin. Phys. 23 077201
[4] Zhang Y R Zhang B Li Z J Deng X C Liu X L 2009 Chin. Phys. 18 3995
[5] Afanasev V V Bassler M Pensl G Schulz M 1997 Phys. Status Solidi 162 321
[6] Lipkin L A Das M K Palmour J W 2002 Mater. Sci. Forum. 389�?93 985
[7] Rudenko T E Osiyuk I N Tyagulski I P Ólafsson H Ö Sveinbjörnsson E Ö 2005 Solid-State Electron. 49 545
[8] Okamoto D Yano H Hatayama T Fuyuki T 2010 Appl. Phys. Lett. 96 203508
[9] Tang X Y Song Q W Zhang Y M Zhang Y M Jia R X Lu H L Wang Y H 2012 Chin. Phys. 21 087701
[10] Tian L X Zhang F Shen Z W Yan G G Liu X F Zhao W S Wang L Sun G S Zeng Y P 2016 Chin. Phys. 25 128104
[11] Iwasaki Y Yano H Hatayama T Uraoka Y Fuyuki T 2010 Appl. Phys. Express 3 026201
[12] Cheong K Y Moon J Kim H J Bahng W Kim N K 2010 Thin Solid Films 518 3255
[13] Fiorenza P Giannazzo F Vivona M La Magna A Roccaforte F 2013 Appl. Phys. Lett. 103 153508
[14] Allerstam F Olafsson H O Gudjonsson G Dochev D Sveinbjornsson E O Rodle T Jos R 2007 J. Appl. Phys. 101 124502
[15] Yoshioka H Nakamura T Kimoto T 2012 J. Appl. Phys. 111 014502
[16] Yoshioka H Nakamura T Kimoto T 2012 J. Appl. Phys. 112 024520
[17] Yano H Hirao T Kimoto T Matsunami H Asano K Sugawara Y 1999 IEEE Electron. Dev. Lett. 20 611
[18] Ueoka Y Shingu K Yano H Hatayama T Fuyuki T 2012 Jpn. J. Appl. Phys. 51 110201
[19] Kimoto T Kanzaki Y Noborio M Kawano H Matsunami H 2005 Jpn. J. Appl. Phys. 44 1213
[20] Nakano Y Mukai T Nakamura R Nakamura T Kamisawa A 2009 Jpn. J. Appl. Phys. 48 04C100
[21] Harada S Ito S Kato M Takatsuka A Kojima K Fukuda K Okumura H 2010 Mater. Sci. Forum. 645�?48 999
[22] Agarwal A K Seshadri S Rowland L B 1997 IEEE Electron. Dev. Lett. 18 592
[23] Singh R Hefner A R 2004 Solid-State Electron. 48 1717
[24] Terman L M 1962 Solid-State Electron. 5 285
[25] Li H F Dimitrijev S Sweatman D Harrison H 2000 J. Electron. Mater. 29 1027
[26] Lee K Y Chang Y H Huang Y H Wu S D Chung C Y Huang C F Lee T C 2013 Appl. Surf. Sci. 282 126
[27] Afanas’ev V V Stesmans A Bassler M Pensl G Schulz M J 2000 Appl. Phys. Lett. 76 336
[28] Rozen J Dhar S Zvanut M Williams J Feldman L 2009 J. Appl. Phys. 105 124506
[29] Rozen J Dhar S Dixit S K Afanas’ev V V Roberts F O Dang H L Wang S Pantelides S T Williams J R Feldman L C 2008 J. Appl. Phys. 103 124513
[30] Nanen Y Kato M Suda J Kimoto T 2013 IEEE Trans. Electron. Dev. 60 1260
[31] Saks N Agarwal A 2000 Appl. Phys. Lett. 77 3281
[32] Šimonka V Hössinger A Weinbub J Selberherr S 2016 International Conference on Simulation of Semiconductor Processes and Devices September 6-8, 2016 Nuremberg, Germany 233 10.1109/SISPAD.2016.7605190
[33] Gupta S K Azam A Akhtar J 2011 Pramana - J. Phys. 76 165
[34] Chanana R McDonald K Di Ventra M Pantelides S Feldman L Chung G Tin C Williams J Weller R 2000 Appl. Phys. Lett. 77 2560
[35] Afanas’ev V Bassler M Pensl G Schulz M Stein von Kamienski E 1996 J. Appl. Phys. 79 3108
[36] Cheong K Y Moon J H Kim H J Bahng W Kim N K 2008 J. Appl. Phys. 103 084113
[37] Arslan E Bütün S Ozbay E 2009 Appl. Phys. Lett. 94 142106